Dynamical glass transition: critical temperature T c and memory kernel in 

MD-simulated Ni 8 Zr 2 
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We use molecular dynamics computer simulations to investigate a critical temperature T c for a dy- 
namical glass transition as proposed by the mode-coupling theory (MCT) of dense liquids in a glass 
forming Nio.sZro.2-system. The critical temperature T c are analyzed from different quantities and 
checked the consistency of the estimated values, i.e. from (i) the non-vanishing nonergodicity pa- 
rameters as asymptotic solutions of the MCT equations in the arrested state, (ii) the g m -parameters 
describing the approach of the melt towards the arrested state on the ergodic side, (iii) the diffusion 
coefficients in the melt. The resulting T c values are found to agree within about 10 %. 



I. INTRODUCTION 

The transition from a liquid to an amorphous solid that 
sometimes occurs upon cooling remains one of the largely 
unresolved problems of statistical physics [H, [3| . At the 
experimental level, the so-called glass transition is gener- 
ally associated with a sharp increase in the characteristic 
relaxation times of the system, and a concomitant depar- 
ture of laboratory measurements from equilibrium. At 
the theoretical level, it has been proposed that the tran- 
sition from a liquid to a glassy state is triggered by an 
underlying thermodynamic (equilibrium) transition Q; 
in that view, an "ideal" glass transition is believed to oc- 
cur at the so-called Kauzmann temperature, Tjc. At Tjc, 
it is proposed that only one minimum-energy basin of 
attraction is accessible to the system. One of the first ar- 
guments of this type is due to Gibbs and diMarzio [J] , but 
more recent studies using replica methods have yielded 
evidence in support of such a transition in Lennard-Jones 
glass formers 0, @. These observations have been 
called into question by experimental data and recent re- 
sults of simulations of polydisperse hard-core disks, which 
have failed to detect any evidence of a thermodynamic 
transition up to extremely high packing fractions [7[ . One 
of the questions that arises is therefore whether the dis- 
crepancies between the reported simulated behavior of 
hard-disk and soft-sphere systems is due to fundamental 
differences in the models, or whether they are a conse- 
quence of inappropriate sampling at low temperatures 
and high densities. 

Different, alternative theoretical considerations have 
attempted to establish a connection between glass tran- 
sition phenomena and the rapid increase in relaxation 
times that arises in the vicinity of a theoretical crit- 
ical temperature (the so-called "mode-coupling" tem- 
perature, Tmct), thereby giving rise to a "kinetic" or 
"dynamic" transition [8j. In recent years, both view- 
points have received some support from molecular simu- 
lations. Many of these simulations have been conducted 
in the context of models introduced by Stillinger and 
Weber and by Kob and Andersen Q; such models 
have been employed in a number of studies that have 



helped sha pe our current views about the glass transi- 
tion |, MM, M M Cl- 
in its simplest ("idealized") version, firstly analyzed in 
the "schematic" approach by Bengtzelius et al. [16| and 
independently by Leutheusser [17]], the MCT predicts 
a transition from a high temperature liquid ( "ergodic" ) 
state to a low temperature arrested ( "nonergodic" ) state 
at a critical temperature T c . Including transversale cur- 
rents as additional hydrodynamic variables, the full MCT 
shows no longer a sharp transition at T c but all structural 
correlations decay in a final a-process [l8[ . Similar effects 
are expected from inclusion of thermally activated mat- 
ter transport, that means diffusion in the arrested state 

Hi- 

In the full MCT, the remainders of the transition and 
the value of T c have to be evaluated, e.g., from the ap- 
proach of the undercooled melt towards the idealized ar- 
rested state, either by analyzing the time and tempera- 
ture dependence in the /3-regime of the structural fluctu- 
ation dynamics [2ll |22| . |23| or by evaluating the temper- 
ature dependence of the so-called g m -parameter [H, [2^| . 
There are further posibilities to estimates T c , e.g., from 
the temperature dependence of the diffusion coefficients 
or the relaxation time of the final a-decay in the melt, as 
these quantities for T > T c display a critical behaviour 
\T — T" c | ±7 . However, only crude estimates of T c can be 
obtained from these quantities, since near T c the critical 
behaviour is masked by the effects of transversale cur- 
rents and thermally activated matter transport, as men- 
tioned above. 

On the other hand, as emphasized and applied in 
[H, H3, H], the value of T c predicted by the idealized 
MCT can be calculated once the partial structure fac- 
tors of the system and their temperature dependence are 
sufficiently well known. Besides temperature and particle 
concentration, the partial structure factors are the only 
significant quantities which enter the equations of the 
so-called nonergodicity parameters of the system. The 
latter vanish identically for temperatures above T c and 
their calculation thus allows a rather precise determina- 
tion of the critical temperature predicted by the idealized 
theory. 
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At this stage it is tempting to consider how well the 
estimates of T c from different approaches fit together and 
whether the T c estimate from the nonergodicity parame- 
ters of the idealized MCT compares to the values from the 
full MCT. Regarding this, we here investigate a molec- 
ular dynamics (MD) simulation model adapted to the 
glass-forming Nio.sZro.2 transition metal system. The 
Ni^Zri-^-system is well studied by experiments I29l . 13011 
and by MD-simulations [H H H, H, \ML HE H3, M, S 
as it is a rather interesting system whose components are 
important constituents of a number of multi-component 
'massive' metallic glasses. In the present contribution 
we consider, in particular, the x — 0.8 composition and 
concentrate on the determination of T c from evaluating 
and analyzing the nonergodicity parameter, the g m {T)- 
parameter in the ergodic regime, and the diffusion coef- 
ficients. 

Our paper is organized as follows: In section II, we 
present the model and give some details of the compu- 
tations. Section III. gives a brief discussion of some as- 
pects of the mode coupling theory as used here. Results 
of our MD-simulations and their analysis are then pre- 
sented and discussed in Section IV. 



III. THEORY 

A. Nonergodicity parameters 

In this section we provide some basic formulae that 
permit calculation of T c and the nonergodicity parame- 
ters fij(q) for our system. A more detailed presentation 
may be found in Refs. [H, E3, H, HO, SI ■ The central 
object of the MCT are the partial intermediate scatter- 
ing functions which are defined for a binary system by 
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N, Nj 



x(e X p(iq-[r l Q (i)-r£(0)]) 



where 



i = 1,2 



(1) 



(2) 



is a Fourier component of the microscopic density of 
species i. 

The diagonal terms a — (3 are denoted as the incoher- 
ent intermediate scattering function 



II. SIMULATIONS 



1 N, 



(3) 



a=l 



The present simulations are carried out as state-of-the- 
art isothermal- isobaric (N,T,p) calculations. The New- 
tonian equations of N — 648 atoms (518 Ni and 130 
Zr) are numerically integrated by a fifth order predictor- 
corrector algorithm with time step At — 2.5 10~ 15 s in 
a cubic volume with periodic boundary conditions and 
variable box length L. With regard to the electron theo- 
retical description of the interatomic potentials in tran- 
sition metal alloys by Hausleitner and Hafner [Z(|, we 
model the interatomic couplings as in [32j by a volume 
dependent electron-gas term E vo i (V) and pair potentials 
4>{r) adapted to the equilibrium distance, depth, width, 
and zero of the Hausleitner-Hafner potentials for 
Ni .gZr .2 33]. For this model, simulations were started 
through heating a starting configuration up to 2000 K 
which leads to a homogeneous liquid state. The system 
then is cooled continuously to various annealing temper- 
atures with cooling rate —dtT = 1.5 10 12 K/s. After- 
wards the obtained configurations at various annealing 
temperatures (here 1500-600 K) are relaxed by carrying 
out additional isothermal annealing runs. Finally the 
time evolution of these relaxed configurations is mod- 
elled and analyzed. More details of the simulations are 
given in [33j . 



The normalized partial- and incoherent intermediate 
scattering functions are given by 



$l(q,t) = F°(q,t) , 



(4) 
(5) 



where the Sij(q) = Fij(q,t = 0) are the partial static 
structure factors. 

The basic equations of the MCT are the set of nonlin- 
ear matrix integrodifferential equations 



F(q,t) + n 2 ( q )F( q ,t) 



drM(q, t — t )F(g, r ) =0 



(6) 

where F is the 2x2 matrix consisting of the partial inter- 
mediate scattering functions Fij(q,t), and the frequency 
matrix Jl 2 is given by 

[Cl 2 (q)] . . = fkBTfa/mi) E s lk [s-\ q )] k . . (7) 

fc 

S(q) denotes the 2x2 matrix of the partial structure fac- 
tors Sij(q), Xi = Ni/N and m, means the atomic mass of 
the species i. The MCT for the idealized glass transition 
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predicts [l8| that the memory kern M can be expressed 
at long times by 

xVifci(q,k)V,7c'i'(q, q- k) 
xF fcfc /(k,t)F„/(q-k,t) , (8) 

where p = AT/V is the particle density and the vertex 
Via/3 (q, k) is given by 

Vi«(q,k) = £«Ci fe (k) + — — — —6 lk cu{q-'k) (9) 



This iterative procedure, indeed, has two type of solu- 
tions, nontrivial ones with f(q) > and trivial solutions 
f(g) = 0. 

The incoherent nonergodicity parameter f"(q) can be 
evaluated by the following iterative procedure: 



Mt[fj:>%) 



(16) 



As indicated by Eq. (JT6J) , computation of the incoherent 
nonergodicity parameter //(<?) demands that the coher- 
ent nonergodicity parameters are determined in advance. 



and the matrix of the direct correlation function is de- 
fined by 



'■'■>(q) = ^-[s- 1 (q)] tf 

Xi J 



(10) 



The equation of motion for F* (q, i) has a similar form 
as Eq.([5|), but the memory function for the incoherent 
intermediate scattering function is given by 



M/(q,t) = / -p- l -(^){cFU^t) 



(2tt) 3 p V q 
Xi^(q-k,t), 



(11) 



B. g m - parameter 

Beyond the details of the MCT, equations of motion 
like (J6j) can be derived for the correlation functions under 
rather general assumptions within the Lanczos recursion 
scheme |44| resp. the Mori-Zwanzig formalism 45]. The 
approach demands that the time dependence of fluctua- 
tions A, B, ... is governed by a time evolution operator 
like the Liouvillian and that for two fluctuating quan- 
tities a scalar products (B, A) with the meaning of a 
correlation function can be defined. In case of a tagged 
particle, this leads for $|(g, i) to the exact equation 



(cF)i(k,t) = (c lt (q)) 2 F tl (q,t) + 2c rl (q)c l] (q)F l] (q,t) 
+ (c lJ (q)) 2 F JJ (q,t) jjLi . (12) 

In order to characterize the long time behaviour of the 
intermediate scattering function, the nonergodicity pa- 
rameters f(q) are introduced as 



/y(q) = Zimt— oo$ij(q,*) 



(13) 



These parameters are the solution of eqs. (|6l)- (|10l) at long 
times. The meaning of these parameters is the following: 
if /ij(q) = 0, then the system is in a liquid state with 
density fluctuation correlations decaying at long times. 
If /y(q) > 0, the system is in an arrested, nonergodic 
state, where density fluctuation correlations are stable 
for all times. In order to compute (q) , one can use the 
following iterative procedure [28| : 



fC'+%) = 



S(g)-N[f(0,f(0](g)-S( g ) 
Z 

, g- a |S(g)||N[f(0,f(0]( g )|S( g ) 



, (14) 



Z = g 2 + Tr(S(g)-N[f« f«](g)) 
+(Z - 2 |S( 9 )||N[fW,f«]( g )| , 

where the matrix N(g) is given by 



XiksT 



Mij(q) 



(15) 



3J(g,t)/ng + *?(?,*) + / drM?(q,t- r)*J(g,r) = 
Jo 

with memory kernel Mf(q, t) in terms of a continued frac- 
tion. 

Within Mf(q,t) are hidden all the details of the time 
evolution of 4>f (q, t). As proposed and applied in (24l.[25|. 
instead of calculating M®(q,t) from the time evolution 
operator as a continued fraction, it can be evaluated in 
closed forms once 3>|(g,t) is known, e.g., from experi- 
ments or MD-simulations. This can be demonstrated by 
introduction of 



$ c (w) ±i$ s M := lim £{$}(e =F«w) , (18) 

e^O 



with 



/>oo 

£{$}(z) = / dte- zt ®(t) 
Jo 



(19) 



the Laplace transform of Q(t), and 



M?(<jj) c ±iM?(w) a := lim£{Mf}(eTiw) .(20) 

£—►0 



Eq.(IT7l) then leads to 



cv > (21) 



On the time axis, Mf(t) is given by 
2 f 00 

M°(t) = - / dujM°(uj) c cos{ojt) 



(22) 
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Adopting some arguments from the schematic 
MCT, Eq. tfTT)) allows asymptotically finite correlations 
$>f(q,t — > oo) > 0, that means an arrested state, if 
M?(q,t-> oo) remains finite where the relationship holds 

MP(g,^co)($|(g,^co)- 1 -l) = l . (23) 

In order to characterize the undercooled melt and its 
transition into the glassy state, we introduced in [24| the 
function 

G($,M°) := M°(i)(l/$(t) - 1) . (24) 

According to IpB"]). G(<f>,M°) has the property that 

G($,Af°) |^oc=l (25) 

in the arrested, nonergodic state. On the other hand, if 

g m := Max {G($, M°) < t < oo} < 1 , (26) 

there is no arrested solution and the correlations $f (q, t) 
decay to zero for t — > oo, that means, the system is in the 
liquid state. From that we proposed [24[ to use the value 
of g m as a relative measure how much the system has 
approached the arrested state and to use the temperature 
dependence of g m (T) in the liquid state as an indication 
how the system approaches this state. 



IV. RESULTS AND DISCUSSIONS 

A. Partial structure factors and intermediate 
scattering functions 



for distance r up to R c = 3/2L. For numerical evaluation 
of Eq. (|27[) . a Gaussian type damping term is included 

Sij(q) = XiSij + AirpXiXj I r 2 (gij(r) ~ 1) 

o 

x exp(-(r/i?) 2 )dr (29) 

with R = R c /3. 

Figffl- figEl shows the partial structure factors Sij (q) 
versus q for all temperatures investigated. The figure in- 
dicates that the shape of Sij (q) depends weakly on tem- 
perature only and that, in particular, the positions of the 
first maximum and the first minimum in Sij (q) are more 
or less temperature independent. 
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01 23456789 10 

FIG. 1: Partial structure factors Ni-Ni-part at T = 1500 K, 
1400 K, 1300 K, 1200 K, 1100 K, 1000 K, 900 K and 800 K 
(from top to bottom), the curves are vertically shifted by 0.05 
relative to each other. 



First we show the results of our simulations concerning 
the static properties of the system in terms of the partial 
structure factors Sij (q) and partial correlation functions 
9ij(r)- 

To compute the partial structure factors S ji(q ) for a 
binary system we use the following definition [46j 

Sij(~q) = Xidij + pxiXj J (gij(r) - l)e~ zq ' r d~r (2,7) 
where 

V I Ni Nl \ 
9 ^ ) = W^X^ ^ 5{v-\v a {t)-r p {t)\)) 

1 i \a=l p=l,pj=a I 

(28) 

are the partial pair correlation functions. 

The MD simulations yield a periodic repetition of the 
atomic distributions with periodicity length L. Trunca- 
tion of the Fourier integral in Eq. (|27|) leads to an oscilla- 
tory behavior of the partial structure factors at small q. 
In order to reduce the effects of this truncation, we com- 
pute from Eq. (|28p the partial pair correlation functions 
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FIG. 2: Partial structure factors Ni-Zr-part at T = 1500 K, 
T = 1400 K, 1300 K, 1200 K, 1100 K, 1000 K, 900 K and 
800 K (from top to bottom), the curves are vertically shifted 
by 0.05 relative to each other. 

To investigate the dynamical properties of the system, 
we have calculated the incoherent scattering function 
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FIG. 3: Partial structure factors Zr-Zr-part at T = 1500 K, 
T = 1400 K, 1300 K, 1200 K, 1100 K, 1000 K, 900 K and 
800 K (from top to bottom), the curves are vertically shifted 
by 0.05 relative to each other. 



F* (q, t) and the coherent scattering function Fij (q, t) as 
defined in equations (JTJ) and ([3]). 




FIG. 4: Incoherent intermediate scattering function 3>f(q,t) 
Ni-part for q = 24.4 nm" 1 at T = 1500 K, 1400 K, 1300 K, 
1200 K, 1100 K, 1000 K, 950 K, 900 K, and 800 K (from left 
to right). 

FigfJ] and figE] presents the normalized incoherent in- 
termediate scattering functions t) of both species 
evaluated from our MD data for wave vector q n =2im/L 
with n = 9, that means qg = 24.4 nm _1 . From the figure 
we see that $f (g, t) of both species shows at intermediate 
temperatures a structural relaxation in three succesive 
steps as predicted by the idealized schematic MCT [47j . 
The first step is a fast initial decay on the time scale of 
the vibrations of atoms (t < 0.2 ps). This step is charac- 
terized by the MCT only globaly. The second step is the 
/3-relaxation regime. In the early /3-regime the correla- 
tor should decrease according to t) — f CS i{q)+ A/t a 
and in the late /3-relaxation regime, which appears only in 
the melt, according the von Schweidler law f CS i(q) — Bt b . 
Between them a wide plateau is found near the criti- 
cal temperature T c . In the melt, the a-relaxation takes 




FIG. 5: The same as fig |3] but for Zr-part. 



place as the last decay step after the von Schweidler-law. 
It can be described by the Kohlrausch- Williams- Watts 
(KWW) law $?(g,t) = A exp(-(i/r ct ) /3 ) where the re- 
laxation time T a near the glass transition shifts drasti- 
cally to longer times. 

The inverse power-law decay for the early /3-regime 
$ ~ f c + A/t a is not seen in our data. This seems to be 
due to the fact that in our system the power-law decay is 
dressed by the atomic vibrations ([24l |25| and references 
therein). 

According to our MD-results, $J(q,t) decays to zero 
for longer times at all temperatures investigated. This 
is in agreement with the full MCT. Including transver- 
sal currents as additional hydro dynamic variables, the 
full MCT comes to the conclusion that all structural 
correlations decay in the final a-process, independent of 
temperature. Similar effects are expected from inclusion 
of thermally activated matter transport, that means dif- 
fusion in the arrested state. 

At T = 900 K - 700 K, the t) drop rather sharply 
at large t. This reflects aging effects which take place, if 
a system is in a transient, non-steady state |48| . Such a 
behaviour indicates relaxations of the system on the time 
scale of the 'measuring time' of the correlations. 



B. Nonergodicity parameters 

The nonergodicity parameters are defined by Eq. l|13|) 
as a non-vanishing asymptotic solution of the MCT 
Eq.®. Fig. [5] presents the estimated g-dependent non- 
ergodicity parameters from the coherent and incoherent 
scattering functions of Ni and Zr at T=1005 K. 

In order to compute the nonergodicity parameters 
fij(q) analytically, we followed for our binary system 
the self-consistent method as formulated by Nauroth and 
Kob [28| and as sketched in Section III. A. Input data for 
our iterative determination of fij(q) = Fij(q, oo) are the 
temperature dependent partial structure factors Sij (q) 
from the previous subsection. The iteration is started 
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Ergodizitaet-Parameter (T =1005 K) 




2 4 6 8 10 12 14 16 2 4 6 8 10 12 14 16 



FIG. 6: Non-ergodicity parameter f ci j for the incoheremt and 
coherent intermediate scattering functions as solutions of eqs. 
© and © 



by arbitrarily setting F Ni - Ni (q, oo) (0) = 0.5Sm-m(q), 
F Zr - Zr (q,oc)W = 0.5S Zr - Zr (q), F Ni - Zr (q,<x>)W = 0. 

For T > 1100 K we always obtain the trivial solu- 
tion fij(q) = while at T = 1000 K and below we get 
stable non- vanishing fij(q) > 0. The stability of the 
non-vanishing solutions was tested for more than 3000 
iteration steps. From this results we expect that T c for 
our system lies between 1000 and 1100 K. To estimate 
T c more precisely, we interpolated Sij(q) from our MD 
data for temperatures between 1000 and 1100 K by use 
of the algorithm of Press et.al. (49|. We observe that at 
T = 1005 K a non-trivial solution of fij(q) can be found, 
but not at T = 1010 K and above. It means that the 
critical temperature T c for our system is around 1005 K. 
The non-trivial solutions (q) for this temperature shall 
be denoted the critical nonergodicty parameters fdj{q)- 
They are included in Fig. [5] 

By use of the critical nonergodicity parameters fdj (q) , 
the computational procedure was run to determine the 
critical nonergodicity parameters f^(q) for the incoher- 
ent scattering functions at T = 1005 K . Fig. [5] also 
presents our results for the so calculated f^(q). 



C. g(3>f , Mj°)-function and g m -parameters 

Here we pres ent our results about the g($f , M- 5 )- 
function [24l. l25j| described in section III.B. The memory 
functions M®(q,t) are evaluated from the MD data for 
t) by Fourier transformation along the positive time 
axis. For completeness, also T = 700 and 800 K data are 
included where the corresponding <l>f(g,t) are extrapo- 
lated to longer times by use of an KWW approximation. 

Fig. [7]and Fig. [5] show the thus deduced M°(q,t) for 
q — 24.4 nm^ 1 . Regarding their qualitative features, the 
obtained M$(q,t) are in full agreement with the results 
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FIG. 7: Ni-Part: Time dependence of the dimensionless 
memory function M° (q, t)/f2s_j from MD simulations for 
q 9 = 21.6 nm" 1 and T = 800 K, 900 K, 950 K, 1000 K, 
1100 K, 1200 K, 1300 K, 1400 K, and 1500 K (from top to 
bottom) 




FIG. 8: The same as fig[7]but for Zr-part. 



in [25j for the Nio.sZro.s system. A particular interesting 
detail is the fact that there exists a minimum in M?(g, t) 
for both species, Ni and Zr, at all investigated tempera- 
tures around a time of 0.1 ps. Below this time, <I?f (<?,£) 
reflects the vibrational dynamics of the atoms. Above 
this value, the escape from the local cages takes place 
in the melt and the /3-regime dynamics are developed. 
Apparently, the minimum is related to this crossover. 

In Fig. |S and Fig. [TO] we display M?(q, t)), 
that means M®(q,t) versus <&| (g, i) . In this figure we 
again find the features already described for Nio.sZro.s 
in 0, [2f| . According to the plot, there exist (q- 
dependent) limiting values <5>f Q (q,t) so that M°(q,t) for 
$f(q,£) < ®to(q,t) is close to an universal behavior, 
while for <5f(<7, t) > $>? Q (q,t) marked deviations are seen. 
$f (q,i) significantly decreases with increasing tempera- 
ture. It is tempting to identify Mf(q,t) below <J>f (g, f) 
with the polynomial form for Mf(q, t) assumed in the 
schematic version of the MCT [18j . In fig. [5] and fig. [TUl 
the polynomial obtained by fitting the 1000 K data be- 
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FIG. 9: Ni-Part: Dimensionless memory function 
M° (q, t)/fig_i as a function of &f(q,t) (solid line) for qg — 
24.4 nm -1 and T = 800 K, 900 K, 950 K, 1000 K, 1100 K, 
1200 K, 1300 K, 1400 K, and 1500 K (from top to bottom); 
a) Ni-part and b) Zr-part ; Polynom fit of the low $ memory 
function M($) at T = 800 K (long dashed line). 
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FIG. 10: The same as fig[9]but for Zr-part. 



low &f (q, t) is included by a dashed line, extrapolating 
it over the whole <E>-range. 

By use of the calculated memory functions, we can 
evaluate the g($f,M°), Eq.p3]). In Fig[ll]and Fig. QJ 
this quantity is presented versus the corresponding value 
of $f (q, t) and denoted as g( < i ) |). For all the investigated 
temperatures, g( < I ) f) has a maximum g m (q,T) at an in- 
termediate value of $. In the high temperature regime, 
the values of g m (g, T) move with decreasing temperature 
towards the limiting value 1. This is, in particular, visible 
in Fig. 1 131 where we present g m (q, T) as function of tem- 
perature for both species, Ni and Zr, and wave-vectors 
qg = 24.4 nm" 1 . At temperatures above 1000 K, the 
g m -values increase approximately linear towards 1 with 
decreasing temperatures. Below 1000 K, they remain 
close below the limiting value of 1, a behavior denoted 
in [24], [25[ as a balancing on the borderline between the 
arrested and the non-arrested state due to thermally in- 
duced matter transport by diffusion in the arrested state 



at the present high temperatures. 
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FIG. 11: Ni-part: MD simulation results for the characteristic 
function g($ s ) as a function of <&1 for q — 24.4 nm -1 
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FIG. 12: The same as figHUbut for Zr-part. 



Linear fit of the g TO -values for Ni above 950 K and 
for Zr above 1000 K predicts a crossover temperature T* 
from liquid (g m < 1) to the quasi- arrested (g m = 1) 
behavior around 970 K from the Ni data and around 
1020 K from the Zr data. We here identify this crossover 
temperature with the value of T c as visible in the ergodic, 
liquid regime and estimate it by the mean value from the 
Ni- and Zr-subsystems, that means by T c = 1000 K. 

While in 0, SH for the Nio. 5 Zr .5 melt a T c -value 
of 1120 K was estimated from g m (T), the value for the 
present composition is lower by about 120 K. A signifi- 
cant composition dependence of T c is expected according 
to the results of MD simulation for the closely related 
Co^Zri-a, system [15j . Over the whole arrange, T c was 
found to vary between 1170 and 650 K in Co x Zti_ x , with 
T c (x = 0.2) ~ 800 K. Regarding this, the present data for 
the Ni^Zri^j, system reflect a rather weak T c variation. 



' 800 1000 1200 1400 1600 1800 

FIG. 13: MD simulation results of the temperature depen- 
dence of g m (q,T) for qg — 24.4 nm _1 (symbols). Linear 
fits to the g m (q, T) are included by full and dash lines ( for 
q 9 = 24.4 nm^ 1 ); a) Zr-part with T c = 970 K and b) Ni-part 
with T c — 950 K. 



D. Diffusion-coefficients 

From the simulated atomic motions in the computer 
experiments, the diffusion coefficients of the Ni and Zr 
species can be determined as the slope of the atomic mean 
square displacements in the asymptotic long-time limit 

(1/Ni) £ K(t)~r a (0)\ 2 
A(T)=lim ^— . (30) 

Fig. HU shows the thus calculated diffusion coefficients 
of our Nio.sZro.2 model for the temperature range be- 
tween 600 and 2000 K. At temperatures above approxi- 
mately 1000 K, the diffusion coefficients for both species 
run parallel to each other in the Arrhenius plot, indicat- 
ing a fixed ratio D^i/Dzr ~ 2.5 in this temperature 
regime. At lower temperatures, the Ni atoms have a 
lower mobility than the Zr atoms, yielding around 800 
K a value of about 10 for DMi/Dzr- That means, here 
the Zr atoms carry out a rather rapid motion within a 
relative immobile Ni matrix. 

According to the MCT, above T c the diffusion coeffi- 
cients follow a critical power law 

D t (T) ~ (T - T c y, for T > T c (31) 

with non- universal exponent 7 [5fj| . 

In order to estimate T c from this relationship, we have 
adapted the critical power law by a least mean squares 
fit to the simulated diffusion data for 1000 K and above. 
According to this fit, the system has a critical tempera- 
ture of about 850-900 K. 

Similar results for the temperature dependence of the 
diffusion coefficients have been found in MD simula- 
tions for other metallic glass forming systems, e.g. , for 
Ni . 5 Zr .5 Ell, for M^Zr^ Cu .3 3 Zro.67 \jm, or 

Nio.8rBo.19 [52(. In all cases, like here, a break is observed 
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FIG. 14: Diffusion coefficients D t as a function of 1000/T. 
Symbols are MD results for Zr (squares) and Ni (diamonds) 

in the Arrhenius slope. In the mentioned Zr-systems, 
this break is related to a change of the atomic dynamics 
around T c whereas for Ni .8iB lg system it is ascribed to 
To- As in 52] T c and Tq apparently fall together, there 
is no serious conflict between the obervations. 



V. CONCLUSION 

The present contribution reports results from MD sim- 
ulations of a Nio.sZro.2 computer model. The model is 
based on the electron theoretical description of the inter- 
atomic potentials for transition metal alloys by Hausleit- 
ner and Hafner (4fj| . There are no parameters in the 
model adapted to the experiments. 

There is close agreement between the T c values esti- 
mated from the dynamics in the undercooled melt when 
approaching T c from the high temperature side. The val- 
ues are T c ss 950 — 1020 K from the g m -parameters, and 
T c « 950 K from the diffusion coefficients. As discussed 
in [1 51 ] . the T c -estimates from the diffusion coefficients 
seem to depend on the upper limit of the temperature 
region taken into account in the fit procedure, where an 
increase in the upper limit increases the estimated T c . 
Accordingly, there is evidence that the present value of 
950 K may underestimate the true T c by about 10 to 50 
K, as it based on an upper limit of 2000 K only. Taking 
this into account, the present estimates from the melt 
seem to lead to a T c value around 1000 K. 

The T c from the nonergodicity parameters describe the 
approach of the system towards T c from the low temper- 
ature side. They predict a T c value of 1005 K. This value 
is clearly outside the range of our T c estimates from the 
high temperature, ergodic melt. We consider this as a sig- 
nificant deviation which, however, is much smaller than 
the factor of two found in the modelling of a Lennard- 
Jones system (28[. The here observed deviation between 
the T c estimates from the ergodic and the so-called non- 
ergodic side reconfirm the finding from the soft spheres 
model[27| of an agreement within some 10 % between the 
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different T c -estimates. 
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